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Abstract: Chiral perturbation theory at finite four-volume V{= L^T) is reconsidered 
with a view towards finding a computational scheme that can deal with any value of M^^L, 
where is a generic Nambu-Goldstone mass. The momentum zero modes that cause the 
usual p-expansion to fail in the chiral limit are treated separately, and partly integrated out 
to all orders. In this way the theory remains infrared finite in the perturbative expansion, 
and the chiral limit can be considered at finite volume. We illustrate the technique by 
computing the quark condensate in a finite volume, smoothly connecting standard results 
in the p-regime for larger masses with those of the e-regime for smaller masses. From 
the partially quenched theory we also obtain the spectral density of the Dirac operator, a 
smooth function from the microscopic region to the bulk region of the p-regime. 
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1. Introduction 

A number of studies have been devoted to the finite volume effects in low-energy QCD, or 
chiral perturbation theory (ChPT) The resulting finite-size scaling theory in ChPT 

has broad theoretical interest as it describes the critical behavior of dynamical symmetry 
breaking, and as such has corresponding applications in statistical physics as well. In lattice 
gauge theory it is certainly of great practical value to have analytical predictions for finite 
volumes available, as they can help in eliminating uncertainties due to the finite sizes used 
in numerical simulations [0. 

To investigate the finite-size behavior of ChPT, essentially two perturbative approaches 
have been proposed so far Q. One is the p-expansion, which has just the same form as 
the perturbative series in an infinite volume, only replacing momentum integrals by the 
discrete sums over momentum due to the quantization in units of 1/L (where L is the linear 
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extent). If we denote by M.^ the mass of a generic (pseudo) Nambu-Goldstone boson, this p- 
expansion is vahd when M-j^L S> 1. It is well known what happens when one takes the chiral 
limit in a volume such that M^L crosses unity and gets even smaller [||, ^: The propagators 
of the pseudo Nambu-Goldstone bosons blow up for one single momentum mode, the one 
of zero four-momentum. This invalidates the usual perturbative expansion, and a different 
technique is required. A solution to this problem was given in |^] in terms of a so-called 
e-expansion. In this scheme the zero-momentum mode is, in a sense that becomes more 
clear below, integrated out exactly, while all the remaining momentum modes are treated 
perturbatively. Since the chiral Lagrangian involves an infinite series of terms, and since 
it is only the perturbative expansion that is jeopardized, "exact" integration here refers to 
the term that is leading order in the quark masses m. 

The e-expansion is thus perfectly suited for studying the extreme case where the quark 
mass m is so small that the pion Compton wave length overcomes the size of the volume, 
MttL <C 1. Since the zero-mode becomes dominant in this e-regime, physical observables 
are mostly dependent on the leading low-energy constants: the infinite-volume chiral con- 
densate S and the pion decay constant F, both in the chiral limit. The next-to-leading 
order terms (with coefficients Lj's) at infinite volume are treated in perturbative fashion. 
This is similar to the p-regime expansion, but the ordering of terms is different. Also in 
this respect, the e-regime provides an intriguing alternative to more conventional ChPT 
since different parts of the chiral theory are being probed to any given order. We note that 
the studies have now also been extended to Wilson ChPT where one has more terms which 
explicitly break the symmetry Q. 

As the e-regime deals with the extremely chiral limit where non-trivial finite-size scaling 
starts to appear (but still far from the symmetric phase since one keeps L ^ 1/Aqcd)) 
universality is at work. Perhaps the most important example of this is the equivalence of 
the zero- mode or vacuum part of the theory to chiral Random Matrix Theory (ChRMT) 
||8|, ^. Little is known in detail on how these universal phenomena cease and the ones 
depending on the dynamics specific to QCD appear when the quark mass increases and 
M-n- becomes of order 1/L [10|. In particular, it is not known precisely how the spectral 
density of the Dirac operator, described by ChRMT in the low end matches on to the 



spectrum at larger scales, in the p-regime [|ri|] -p3||. 

Recent developments in both computational facilities and algorithms have allowed 
simulations of full lattice QCD near the chiral limit, but no study has until now reached 
deep inside the e-regime except at rather strong coupling. Although results have often 
compared favorably to the e-expansion of ChPT, there may still be large systematic errors 
due to the condition M.,^L <C 1 not being well fulfilled (see, e.g., ref. [0]). One might 
therefore ask whether it is possible to have a new approach which smoothly connects the p- 
expansion and e-expansion and which remains valid even in the region M^L ~ 1. Recently, 
steps have been taken in that direction by means of a so-called mixed expansion |15, 16| 
(see also ref. ||T7[| ), where one treats the very light flavors with the counting rules of the 
e-expansion, while heavier flavors are counted according to the p-expansion. The results 
turned out to be mixtures of the properties of the e and p regimes: zero-mode fluctuations 
from the light sector in addition to 1-loop corrections from the heavier sector that include 
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chiral logs and some of the -Lj's. But Refs. [15, 16] treated the light and heavy flavors 
separately and did not attack directly the regime where M^L ~ 1. Actually, two regimes at 
play here: one is the first obvious threshold when M^^ ~ the other is when M^^ ~ l/L?'-, 
the scale of the e-regime. The question is what happens in-between. 

In this paper, we suggest a new perturbative approach where all the terms in the p- 
expansion are kept but the zero mode is treated in exactly same way as in the e-expansion^ . 
The expansion thus considers the zero momentum mode on a different footing from the 
rest, partially resumming terms to all orders. Before reaching the e-regime this means that 
an infinite series of terms that are normally considered in the perturbative expansion are 
included to all orders. The result is a slightly re-ordered perturbation theory expansion 
that is free from perturbative infrared singularities in the chiral limit. As an example, we 
compute here a formula for the chiral condensate which smoothly connects the results of 
the p-regime PI and the e-regime. We will argue that our formula is reliable even in the 
intermediate region where we go from M^^ ~ 1/L to Mj^ ~ l/L?. 

Using the partial quenching technique based on replicas, we can treat a general theory 
with valence flavors and physical sea quark flavors of masses which are non-degenerate. 
This allows us to take the discontinuity on the imaginary axis of the valence quarks, and 
thus obtain the corresponding spectral density of the Dirac operator. It is also given by a 
smooth formula that connects known results in the e-regime and p-regime [0, 12 1. 

The rest of our paper is organized as follows. In Section ^, we describe in detail our new 
perturbative method in ChPT. The chiral condensate to next-leading order is calculated in 
Section ^. Taking the discontinuity on the imaginary axis, we obtain the spectral density 
of the Dirac operator in Section ^. We show in Section ^ that our results provide a 
smooth connection between the e and p regimes. In Section |6| we present some numerical 
examples which are useful when comparing with lattice QCD simulations. Conclusions and 
an outlook are given in Section |^. 



2. A chiral expansion at finite volume 

Let us consider A'^^-flavor chiral perturbation theory in a finite volume iV = L^T), 

C = ^Ti[d^U{x)'^d^,U{x)] - ^TT[M'^e-'^I^W{x) + U{x)^ e'^'^f M] + • • • , (2.1) 

where U{x) G SU{Nf) and 6 denotes the vacuum angle. Here, T, is the chiral condensate 
and F denotes the pion decay constant both in the chiral limit. There are of course next- 
to- leading order terms, indicated here by ellipses, with additional low-energy constants 
denoted by Lj's, ffj's and beyond. 

In the partially quenched case, we use the replica method where the calculations are 
done within an (Nf + + (N — A'^y))-flavor theory followed by the replica limit — > 
H^, The ordinary physical Nf -theory result can clearly be viewed either as one where 

^We understand that F. Niedermayer (unpublished) has considered an analogous scheme in the context 
of the 0(n) sigma model. 
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N = = 0, or, alternatively, one where m„ = mf with nif denoting one of the physical 
quark masses. 

From now on we consider sectors of fixed topology u, obtained by Fourier transforming 
in 9 in the usual way. This extends our integration from SU{N + Nf) to U{N + Nf) in 
the zero- momentum sector. 

For the mass matrix, we consider a general diagonal case, 

M = diag(m„, mt,, • • •, mi, m2, • • ■), (2.2) 

N Nf 

where we have = + {N — Ny) replicated flavors and A'^^ physical flavors. Here 
A'^ is the number of, in this case, degenerate valence quarks. What we do below can 
straightforwardly be generalized to non-degenerate valence quarks by just adding copies of 
each. 

We start by factorizing the flelds into the zero-momentum mode U and non-zero modes 

U{x) = C/exp(iV2C(x)/F), (2.3) 

and expand perturbatively in S^{x) just as in the e-regime Q. But here we give the same 
counting rules for the fields and other parameters as in the p-regime: 

a^~0(p), e(x)~0(p), M^O{p^), T,L^O{l/p), (2.4) 

in units of the cut-off AttF. The aim is to see if we can tune quark masses so that we go 
from mass scales M^O{l/L^) through 7W ~ ©(l/L^) to zero. Here we of course assume 
that the linear sizes of the volume, L and T, are much larger than the inverse QCD scale 
^QCD so that the effective theory is valid. 

The above parametrization Eq. (2.3) leads to a well-known Jacobian in the functional 
integral measure [^. Although it is easily taken into account, its contribution is 0{p^) and 
beyond the accuracy with which we do actual calculations in this paper. 

We now expand the Lagrangian in ^{x) according to the p-counting Eq. ( |2.4D , and 
write down the terms relevant to one-loop order for the chiral condensate. 



C = --Tr 



+^Tt[MHu - l)e + fiU^ - 1)M] 



2 



-F2(p) Tr[A^t^] + ..., 



(2.5) 
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where Mfj = (mj + m)T,/F and Lj's are the usual higher-order low-energy constants of 
ChPT. Here we have added and subtracted a conventional mass term of the p-regime. 



We will treat the mass term of the first line in eq. (2.5) as part of the exact Gaussian 
integration that leads to the conventional massive propagator of the p-regime, while the 
remaining terms are treated in a perturbative expansion. We return to this point below. 
The contact term H2 has no direct physical significance, but it is needed as a counterterm 
for the one-loop correction to the condensate All linear terms in ^ are absent due to 

J d^x Cix) = 0. (2.6) 

In Eq. ( |2.5| ) the first line contains terms that in the usual p-expansion are of order 
(the first, a trivial constant in the usual infinite-volume p-expansion) and (the remaining 
two). In the e-expansion the first two terms on the same line are of order e^, while the 
third is of order e^. In the last four lines we have written out explicitly those terms that 
are of order in the usual infinite- volume p-expansion (but trivial constants there) . In the 
e-counting these terms are of order e^. Other terms involving the Lj's will be of order p^ in 
that same counting. However, as with the measure term, these terms will not contribute 
to the chiral condensate that we will compute below. 

The one single term we have not yet discussed is that of 

A£ = ^Tr[M\U - l)e + fiU^ - 1)M] . (2.7) 

In the usual p-expansion a term of this type first occurs at order p^ (because there will 
be three powers of ^), and in the e-expansion it is of order e^. Here we treat any matrix 
elements of Ai^{U — 1) (and its complex conjugate {U^ — as of 0{p^) for all values of 

Ai. By performing the exact group integration over U, we can check that the combination 
A4^(C/ — 1) gives NLO contributions (< 0{p^)) to the results. We illustrate this in Section 
|5[ In Appendix^ we describe an alternative method which expresses the magnitude of the 



contribution from the term in eq. (2^) directly in terms of masses and the volume V, thus 
giving a precise counting of this term in terms of p ~ 1/L. This alternative method gives 
identical results, but is in practice more cumbersome than the scheme presented here. 
By taking 7Wt(C/ - 1) ~ 0{p^), one can thus, to this order, rewrite the Lagrangian 



-Tr 
2 



+^Tt[MHu - l)e + fiU^ - 1)M] 



-L, (^§Tr[M^U + U^M]^ ' - ^^^^ (p) ' Tr[{M^U + U^Mf], (2-8) 

where the second line is treated as an NLO interaction term and the contribution from the 
Ly-term has been dropped. 

It should be stressed at this point that adding and subtracting an ordinary p-regime 
mass term and then expanding the term ( |2.7D perturbatively has the effect or a complete 
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re-ordering and partial resummation of the perturbative series. This resummation comes 
from the fact that when the chiral hmit is taken and is no longer of order p^, but 
smaller, we still keep the full massive propagator. The error in doing this, rather than 
expanding the propagator to the needed order in M^, is however always of yet higher order 
and part of the unavoidable uncertainty in any fixed-order perturbative calculation. We 
always keep the full massive propagator in the expressions and plots we present below. 

The Feynman rule for the .^-propagator is thus obtained as usual, except that the 
zero-momentum modes are not included: 

{^i3{x)^ki{y))!i = Su6jkA{^ - y, Mf.) - 6ij5MG{x - y, M/., m|J, (2.9) 

and the second term comes from the constraint Tr^ = 0. The propagators A and G are 
given by^, 

1 e^P^ 
G{x,MlMf^) = ^- ^, (2.11) 

P^O {p2 + M2. ) (p2 + M2 . ) E V ' 



where the summation is taken over the non-zero 4-momentum 

p = 2TT{nt/T,nx/L,ny/L,nz/L), (2.12) 

with integer n^. 

3. The chiral condensate 

The chiral condensate of a valence flavor is obtained in the conventional manner by adding 
a source to the mass matrix; M ^ A4 + and differentiating the partition function 
Z{M + J) with respect to J . To leading order in our expansion this gives 



(9.9.)^^ ^ i^ln^(A^ + J) 



= ^(f/.. + C/i)^ , (3.1) 
^=0 ^ ' ^ 



where the zero-mode integral (iJvv + ^vv^ is computed non-perturbatively with respect 
to the zero-mode partition function 

dU (dett/j-'exp 

/(7(Af+Af/) 



Zlo = [ dU (detZ7)''exp 



(3.2) 



The analytical formula is known for the most general partially quenched case with non- 
degenerate physical A^j-flavors [ 21 1 . Some details are summarized in Appendix |^. Here we 
simply define 

tl'^{|l,,{^isea]) ^ \{Uvv + Ul)^ , (3.3) 



■^We do not consider the fully quenched theory in this paper. We thus have Nf 7^ in all that follows. 
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where ji^ = m„Sy and the set of the dynamical flavors are denoted by {/J^sea} = tJ'2, • • •} 
with Hi = niiTiV. 

At next-to-leading order, it is convenient to first calculate the 1-loop perturbative 
correction due to the non-zero modes. This can be done by simply evaluating 



1 - J d^x^TriiM + J)\U - l)e + eiU^ -l)iM + J)] 



(3.4) 



where {■■■)£, denotes the integral over ^, and then re-exponentiating it. The effective 
Lagrangian (with a scalar source then reads 



(3.5) 



where 
Zi = 1 



1 

F2 



A(0, M2 ) - G(0, Ml M|) - 16L6 Mf^ - A{2Ls + F2)M| . (3.6) 

j j 

Both of the 1-loop integrations A(0, M^) and (5(0, M^, M^) are UV divergent and their 
divergences are absorbed into the bare parameters L^, Lg and H2. With an appropriate 
regularization such as dimensional regularization (see ref. Q for a discussion of this issue) 
A(0,M2) is given by 

M2 



A(0,M^ 



167r2 



(lnM^ + ci) + (7i(M^), 



(3.7) 



where ci represents the conventional logarithmic divergence which is independent of M 
and the volume 0]. The function gi represents the finite size effects. 



51 



1 



(4vr) 



Edr 



-rM^ 



-yj 



1 



M21/' 



(3.^ 



and the sum is taken over a 4-dimensional vector = n^L^ (L^ being the lattice size in 
the H-th direction) with integer n^. It is particularly important to note that the massless 
limit — > of ^1 is finite and given by 51 (0) = —/3i/VV where Pi is the so-called 
shape coefficient [^. This term dominates the NLO correction in the e-regime. A detailed 
numerical treatment of A(0, M^) and its derivative is discussed in Appendix |C|. 
Expanding G{0,M'^,M'^) in (see Appendix |D|) , we get 

G(0, M2, M2) = A A(0, M^)-^Y1 ^(0' ^JJ^ + • • • ' (3.9) 



N 



f 



f 



where the sum is taken over physical flavors only. Since the O(M^) terms indicated 
by ellipses are UV finite, the divergence in Zi, Eq. (|3.6D, is in total, 



Z,. 



1 



167r2F2 



2 

j 



Cl, 



(3.10) 
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which can be absorbed into redefinitions of Lg, Lg and H2, 




H2 = ( ^ - — ^ ci + m, (3.11) 

(167r)2 \ 2 Nf J ^ ^ ^ 

where the renormahzed constants are denoted by Lg, Lg and i?2- This renormahzation is 
identical to that of the infinite volume case Hill, as it should be. 



To this order, the chiral condensate at fixed topology can thus be written 



(3.12) 



where fi'^ = Z^m^TM and {//^eal = {ZimiT,V, Z2m2'SV, ■ ■ •}. Note that the arguments of 
the function now include the chiral logarithms as explicitly seen in Zi (see Eq. ( |3.6| )). 

The expression in Eq. ( 3.12| ) looks simple and compact. But in order to see the valence 
mass dependence, it is more convenient to decompose Zy into two finite pieces. 



Z„ = Z° + 6Zv{my), 



- 1 ^ 



A(0, Mfj/2) - G{0, 0, 0) - 16L6 M: 



(3.13) 
(3.14) 



V(A(0,M?)-A(0,M,V2)) 



j 



r2 i\^2 



(G(0, M^„M^,) - G(0, 0, 0)) - 4(2L8 + i?2)M, 



^ 1 



where only 5Z^{my) has a valence mass dependence, and it vanishes in the limit m. 
Note that G(0, 0,0) is infra-red finite (see Appendix |d|) . 

With the above decomposition, the condensate can be expressed as 



(3.15) 
, ^0. 



{qvqv)u 



d 



Z0sPQ(Z0/i„,K,J)+5Z,(m,)l +0(/) , 



(3.16) 



where fiv = niv'^V. In the second line we have used the fact that 6Zy{mv) ~ 0{Myy) does 
not contribute until f^u^d^y, {^'g^^}) becomes close to 1 + 0(^2). By the same technique, 
one can replace the sea quark's argument fi[ = ZirriiTiV by Z'^miTM without producing 
no additional term. The explicit form will be given in Section 6. 

It is now clear that the chiral condensate near the chiral limit is dominated by zero 
modes, and hence expressed through combinations of Bessel functions as in the e-regime. 
The argument Z^myTM ^ however, includes the chiral logarithm of the sea quarks in Z^. 
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As the valence mass increases, {/isea}) approaches unity and the ordinary valence 

quark chiral logarithm appears in 5Zy{mv). For yet larger values of m^,^, the term pro- 
portional to 2Lg + H2 becomes important. The unphysical quantity H2 depends on the 
regularization scheme, and the condensate is then not unambiguously defined, as is well 
known |^]. In such a region, one has to eliminate the H2 dependence to obtain unam- 
biguous physical observables. An example would be to consider a difference between two 
topological sectors. As we will see in the next section, the spectral density is also free from 
this ambiguity. 

4. The spectral density of the Dirac operator 

In the previous discussion, we assumed that all quark masses were real and positive. As 
is well-known Q, by considering the expressions for imaginary valence quark masses, one 
can calculate all spectral correlation functions and individual eigenvalues distributions of 
the Dirac operator of the A'^^-flavor theory. In this case partial quenching is simply used to 
extract a physical observable in the full theory. One expresses the valence quark condensate 
as a spectral sum over the Dirac eigenvalues iA^'s (Afc is real), 

= ^y] (--^r-v-) ■ (4-1) 



V ^ \ + iXk 

Since every non-zero eigenvalue comes paired with one of opposite sign, the condensate 
satisfies 

{qvqv)\m* = {{qvqv)\rnj* , {qvqv)\-m, = - {qvqv)\m^ , (4.2) 

where * denotes complex conjugation. 

Using the above, the spectral density at fixed topology v is given by 



Elim 



2vry ^e^o\i(A-FAfc)-e i(A-hAfc)-he 

= lim Tr{'\qvqv)v\m^=i\-f. — '\qvqv)v\mv=%\+^) 
e— »0 ZTT 

= - [Z° ReEPQ(iASyzO, {/.',,J) + Re {bZ^X)^ , (4.3) 

where we have neglected the (5(A) term which represents the exactly zero eigenvalues due 
to the non-trivial topological charge v. Similar expressions exist for all higher spectral 
correlation functions. 

For the calculation of Re (5Zt,(zA)), we need the discontinuities of various functions. 
For the log-terms, for example, one obtains 

Reln(m+ v^)U„=iA = ^In(m2 + A2), (4.4) 

Imln(m -|- \Anf)L„=jA = arctan — , (4.5) 

m 



'Because m„ — /ii,/(EI/*), such a term first appears at NNLO in the pure e-expansion. 
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where m denotes a real and positive mass, for which the hmit m — > is well-defined. Some 
of the other functions occurring in eq. ( |4.3| ) need to be treated numerically, such as those 
in gi. We collect some representations useful for numerical purposes in Appendix Note 
that the term proportional to 2L8 + H2 has disappeared upon taking the discontinuity 
across the imaginary valence quark axis^. 

Next, using the properties of Bessel functions, one sees that the first term of Eq. ( [4.3| ) 
reproduces the known form of the leading contribution to the microscopic spectral density 
in the e-regime, 

ReSPQ(iC, {//sea}) = Vrp--(C, {^^sea}), (4.6) 

where /5™'=(C, {/i^ea}) is given by @ 

Id detB 



2n7(c2+Mf)d«t^ 

Here the Nj x Nf matrix A and the {Nf + 2) x (Nj + 2) matrix B are defined by 



(4.7) 



Aij 


— n^^^T 


-j-lifJ'i) 








= C-'Ju- 


h,-2(C), 








= (-Mi-2 




-l(^i- 


-2) (^/1,2). 



(4. 



(4.9) 

The general formula for the spectral density can thus conveniently be written in a 
representation suitable for small eigenvalues that go into the "bulk" region, 

p,{X) = - [kZI pT^ZIXLV, {^',, J) + Re {bZ,{i\))\ . (4.10) 

TT 

5. Prom the p-regime to the e-regime 

In this section we explain how our formulae are consistent with known results in both e 
and p regimes. We next consider the validity of our expressions in the intermediate region. 

5.1 Checks in the p-regime and limit to the e-regime 

We first check that the expression derived in the previous section reproduces known re- 
sults in the conventional perturbative p-expansion. The exact zero-mode integral above 
is expressed by complicated combinations of Bessel functions. But for large rriiTiV, these 



Bessel functions can be expanded in l/rrii'EV (see, e.g., ref. |23] for details), 

EPQ(/X„ {p,ea}) = 1 - E + + ^^P'^' (^-1) 

j f^v + Pj 8pf, 



*This result may seem to contradict the fact that the condensate, conversely, should follow from a 
spectral sum involving p(A) as in Eq. ([4.l|). The problem is that the spectral sum is UV divergent. It is 
this divergence that gives rise to an ambiguity such as indicated by the 2Ls + H^-term. 
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which, after summing over topology, gives 



+ 



(5.2) 



where we have used 



+ 



(5.3) 



in the chiral limit. Substituting Eq. ( |5.2| ) into Eq. ( |3.12| ), we reproduce the perturbative 
p-regime result, 



{qvqv) 



2 ^ 



(5.4) 



As is known from the matching between e and p regimes, the zero-mode fluctuations which 
give rise to the second and third terms in Eq. (|5.2| ) are absorbed in the momentum sum in 
A and G so that we recover the usual propagators 



G{x,Ml,M. 



J3^ 



-y 



P (p2 + M^)(p2 + M2.)(E 



1 



(5.5) 
(5.6) 



of the ordinary p-expansion. 

The above expression agrees with known results that can be found in, e.g., the work 
of Osborn et al. They derived a formula for the partially quenched case with Nj 

degenerate flavors of mass M"^^^ in which case one can use 



G{x,Ml,Ml) = ^ [A(x,m2J + {Ml - mD5a/2^A(x,m2J] 



(5.7) 



One can also check the e-regime results at fixed topology v are precisely reproduced 
just by reducing the quark masses in the formula Eq. ( 3.16| ). One notes that in that limit. 



1 + 



6Zy{mv) — > 0. 



(5.8) 
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The spectral density of the Dirac operator can easily be compared to known results in 
different limits as well. First of all, the Banks-Casher relation is trivially reproduced when 
we take the limit F — > oo before nii — > 0, 

lim lim p^{X = 0) = -, (5.9) 

rrii^OV^oo TT 

after having used /5Jf*'^(oo, {/Ugea}) = ^/t^, Z^\v^oo,m^^o — > 1, and 6Zy{0) = 0. In the same 
limit above, but for finite A, one obtains 

lim lim p^{X) = — [1 + Re6Z^{iX)\v^oo,mi^o] 



jTii— +oy— >oo vr 
_ S 
vr 




327rF4 ' 



(5.10) 



which is the result of Smilga and Stern |11|. If one keeps the sea quark masses nii = m 
finite (and degenerate), one gets 

lim p^{X) = — [ZqIv^oc +'^e6Zy{iX)\v-,oo] 

V^oo TT 



1 + ^^W^TTT^. ( 2iV/|A| arctan ^ - 4vr|A| 



m 



TT 



(5.11) 



where iXsub denotes the subtraction scale. This is consistent with the formula by Osborn et 
al. [jl2|. In the case of finite V and very small and A, the general result in the e-regime 
11, 

/9,(A) = Spr^(ASy,/x,ea), (5.12) 
is easily recovered upon noting that to leading order we have Z^ = 1. 
5.2 Intermediate regime 

As seen in the above discussion, our formulae for the condensate and the spectral density 
smoothly connect the results in the e-regime with those of the p-regime. But we also need 
to know the precision in the intermediate region. The e-regime assumes M^^L <^ 1 while 
the p-regime counting requires M^^L ^ 1. Our re-ordered perturbative expansion has 
removed this constraint. Instead, we make use of a non-trivial prescription for the fourth 
term in Eq. ( |2.5| ) where we take Ai^{U — 1) to be always small, specifically of 0{p^) or 
smaller. 

One should note that non-zero mode's contributions are free from infra-red divergences 
by construction. This is also seen explicitly in the finite chiral limit of of the finite- volume 
function g{M'^). For the non-zero modes, there is no need to distinguish the e-regime from 
the p-regime. The smaller the quark masses, the better convergence of the non-zero mode 
expansion. 
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Therefore, the accuracy of our calculation needs to be assessed by considering the zero- 
mode integrals. To show the general validity of the method, we need to confirm by explicit 
evaluation of the group integrals that the operator M^{U — 1) consistently can be taken 
to be of 0{p^) or smaller. This in any combination of matrix elements, power and for an 
arbitrary choice of the mass matrix, including such as is needed for partial quenching. 

However, for the calculation of the chiral condensate in this paper, we only need to 
check in Eq. ( p| ) that the U integral keeps the second term of NLO. Since {^fj) = 
it is in fact enough to confirm that [Ai^{U — 1) + h.c.ja ~ 0{p^) or 

m,(^^{Uu + ul)u-l^ -Oip"), (5.13) 

for any i, which can be done directly by means of the exact group integration Eq. (p.3|). 
Without using the rather complicated exact expression, its asymptotic behavior is known 
|23] for niiTiV <^ 1 and niiTiV ^ 1, and this leads to 



rrii 



2 ; 1 _ (AT^ _ _L) _ (rrii^V oo) 



with the other masses mj's (j ^ i) fixed, where nj denotes the degeneracy of the mass rrii 
(Note that n, = in the partially quenched case) . Since this function is everywhere regular 
for finite rriiTiV, one expects that the two limiting cases above are smoothly connected and 
the function thus always kept small, here of 0{1/V) ~ C'(p^) ^■ 

In Fig.Q, we plot the function niiT, (^{Ua + Uj-)u/2 — 1^ for various cases in a (2 + 1)- 
flavor theory. Every curve shows a monotonous function connecting the two limits, thus 
confirming that the Ai{U — 1)^^ contribution to the condensate is always of order 1/V. 
We provide an explicit analytical expression in an analogous U{1) toy model in Appendix 



6. A few examples 

In this section, we present two explicit numerical examples. One is an Nf = 2 degenerate 
two-flavor theory and the other is an Nf = 2 + 1 theory including a strange quark whose 
mass is different from up and down quark masses. For the low-energy constants in both 
cases, we take the following phenomenologically reasonable values: S^/^ = 250 MeV, F = 
90 MeV, Llinsub = 0.77GeV) = 0.05 x lO'^ and LKusub = 0.77GeV) = 0.5 x 10'^ ^j^g^g 
/isub denotes the subtraction scale. 

6.1 degenerate Nf = 2 

Let us first consider the two-flavor theory with degenerate up and down quark masses 
rriu = = m. The factor = + 5Zy{my) is then explicitly given by 

1 



Z?, = l 



%\n^+2-g, {Ml/2) 



^We note that M{U — 1) contributes not as 0{p^) but even further suppressed, of C(p*). This is due 
to the fact that we are here considering a one-point function. 
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1 

+ 2 



Pi 



K'(Y^lnl^^/VL6 + /32)}-32L^(/.,„,)M2 , (6.1) 



and 



Ml 
167r2 



V 2 



m 



1 + ^ In 1 + — + 



m 



m 



167r2 



In- 



Ml 



m 



2/i 



2 

sub 



V 2 



m 



A 



m 



167r2 



52 M,^-f +-.ln 



87r2 



(6.2) 



where M^ = {rriu + mdjTi/F'^ = 2mT, / F'^ .Foi the numerical implementation of gi and g2, 
see Appendix Note that 6Z^{my = 0) = 0. 
For the spectral density, we also need 



Re5Z^{iX) = -^2 



Ml 
167r2 



1 



A2 



A 



A 



m 



-2Re5i 1 + - 

2 V m 



— In ( 1 H ^ I arctan — 

m 

1 

' 2 



2.. I ^ 



Reft ( A/^^^) + * 



A 



+Af^- ( Im52 [Ml-]- — 



iX 



m 



1 



167r 



(6.3) 



As discussed above, the H2 (and /iswfe) dependence has disappeared. 



The non-perturbative expressions for the zero- mode integrals are given by [21| 



det 



(^2_;,2)2 

^fl^{^4Ku+2{^J'v)) fJ-llu+2{fJ-v) fJ''^Iu+2{fJ') fJ'Iu+lifJ') 

\-^^lv{^J'lKy+z{^^v)) lAiu+ziiJ'v) ^J!^Iu+?.{^J) n^iu+2{fj') J 



det 



(6.4) 



and 



ICI 



2(/i2 + C2)2 
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det 



C^z.+l(C) C'^Ju+2{0 lJ''^Iy+2{p) fJ-L+lil^) 



V CV,+2(C) C=^J.+3(C) -/U=^/.+3(m) -^'4+2(m) / 



det 



where /i^, = m^TiV, fi = mTiV and (" = XSV. 



(6.5) 



Substituting the exphcit expressions above into Eq. and (p^ ), it is straightfor- 

ward to compute numerically the chiral condensate and the spectral density of the Dirac 
operator. Using the numerical values listed in the beginning of this section, we plot the 
chiral condensate and the spectral density in Fig. ^ for the case of niu = md = 10 MeV, 
L = T/2 = 2 fm, 1/ = and here taking just as an example i?J(0.77GeV) = 0.1 x lO'^. 
One sees that the resulting thick line smoothly connects the one in the e-regime (solid) 
with the one in the p-regime (dotted)^. It is clear that the effect of one- loop corrections 
to the spectral density of the Dirac operator in the p-expansion is relatively small when 



Nf = 2, in agreement with the prediction Eq. ( ^.10 ) of Smilga and Stern |11] 



It is useful to see what happens when we vary some parameters. In the larger volume 
{L = T/2 = 3 fm) shown in Fig. |3| the agreement between the three different expansions 
becomes visibly better around ML ~2-3. The sea quark mass dependence is shown in 
Fig. § where we show the plots for different sea quark masses at 2 MeV (solid), 10 MeV 
(dotted), and 30 MeV (small dotted). The black filled circle in the left part shows the 
physical points where rriy = ruu = nid- In Fig. |5| the z^-dependence is presented at a fixed 
value of rriu = 10 MeV and L = T /2 = 2fm. The topology dependence disappears around 
MyyL = 3-4. The condensate of the physical theory with m„ = = is plotted on the 
left part of Fig. ^ after summing of topology (See Appendix]^ for the details). Of course, 
the condensate in the massive case is inherently ambiguous due to the presence of the 
coupling H2. On the other hand, the spectral density is free from this ambiguity. On the 
right part of Fig. |^ we plot the density after summing over topology at fixed m„ = = 10 



MeV. As has been observed before |24], the striking oscillations that are seen at sectors of 
fixed topological charge become smeared out upon summing over topology. In the massless 
limit the oscillations will reappear and the density approaches that of the = sector. 



%.2 Nf = 2 + 1 

Next, we consider the A'^^ = 2 + 1 theory where the strange quark mass is different from 

® Because we do not have p-regime predictions available at fixed topological charge u, we always compare 
with p-regime expressions where topological sectors have been summed over. Of course, p-regime predictions 
at fixed topology will differ slightly from these, but the difference becomes insignificant at large volumes. 
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the up and down quark masses. Then G{x, M^, M^) becomes much more comphcated, and 
we refer to the exphcit expressions given in Eqs. ( plsD and (pl6|) . 

One obtains 



1 



A/f2 A/f2 /\/f2 

^ ^ + (M,y.) + ^ , (Mi/2) 



+ 1+ 



16L^(/x,„,)(2M2, + m2 



M 

W) Ml, + 2M2 



+ 91 



1 



(6.6) 



and 



1 

+251 ( 



167r2 



In 



1 + — 

m, 

1 + — 



In 



nij, 
1 + — 



+ 



M2 



2/x2 



+ 51 



-251 



167r2 



' 2/^2 



1 + ^ 



327r2 
-51 



sub 
1 + — 

m 



M2 



' M2 



■ 2A-4 



3 X 167r2 



In 



Ml, + 2M'^ 



In 1 + 



3^ 



■+51 



^+5i(M, 



2 N 



1 + 



327r^ m 
2{Ml, - Mlf 
mi - Ml, - 2Mlf 

mia - Mir , , __ 

^ V (3M2„ - Ml, - 2Mlf) Vl67r2 

ml - Ml,){Ml - Ml) f 1 ..Ml 2 
+ 3M2,-M,^,-2M^ ll^^^^S ^ 1) +5m.^i(M., 
, mid- Ml? ( Ml, + 2Ml Ml, + 2Ml 
^{Ml, + 2Mlf \ 3xl67r2 
2{Ml,-Ml)^ \ ( /5i 
(^L + 2M2J2 



A^d + 



2M2, 



3^^ 



+ 51 



Ml, + 2M2 



-4(2L^(^,„6)+/72'(/^.«6))M, 



JMldMl^f 1 i,yi/2,2 . 



+ 



ot] ' 



(6.7) 



where its real part for the imaginary valence mass is taken as in the same way shown in 
the previous subsection. Here Ml, = {mu + md)Y./F'^,Ml = 2msJ:/F'^,Ml = 2m^T./F'^. 

The non-perturbative zero-mode integrals are given by 
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/ 



det 



dt,AfJ-lKu+2{f^v)) flllu+2{flv) f^llu+2{f^u) HuIu+lifJ'u) fi%+2{fJ's) 
-df^Af^lKi^+sif^v)) ullu+ziHv) fllL+sif^u) /i^4+2(/in) fl^I^+sifJ's) 
df,^{flf,K^+4il^y)) tljl^+iilJ'u) Hjliy+sil^u) H%+A{fJ's) / 



det 



(6. 



where = niu^d'^V and Hs = rrigTiV and one can take the imaginary value of = to 
obtain ^Ug}) in the same way as the 2-flavor case. 

Substituting again the numerical values in the beginning of this section, we plot the 
chiral condensate and the spectral density in Fig. |^ for the case with = = 10 MeV, 
rus = 110 MeV, L = T/2 = 2 fm, u = and iJ^(0.77GeV) = 0.1 x lO'^. Due to the 
large contribution from the chiral logarithm of strange quarks, both of the (normalized) 
condensate and spectral density shows a quite larger deviation from 1, which is the universal 
chiral limit when V ^ oo with any number of flavors. 



7. Conclusions 

We have computed the quark condensate and the spectral density of the Dirac operator by 
means of the chiral Lagrangian in a finite volume. We have used a new perturbative method 
which keeps all the terms which appear in the conventional expansion in the p-regime but 
which treats zero-mode integral non-perturbatively in the same manner as in the e-regime. 
The resulting perturbative series connects smoothly previous results in the p-regime with 
those of the e-regime. Having analytical formulas available in the intermediate region may 
be of obvious value for lattice QCD simulations. 

It would be interesting to investigate the present proposal for finite-volume chiral 
perturbation theory in greater detail. One should obviously consider using such a method 
to compute other finite- volume observables in chiral perturbation theory, such as two-point 
correlation functions. 
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A. Non-perturbative zero- mode integrals 



In this appendix, we briefly review how to perform the non-perturbative zero-mode in- 
tegrals. We are interested in the most general partially quenched calculations for both 
practical purposes of comparisons to lattice data and for computing spectral correlation 
functions of the Dirac operator. Of course, results for the A'^^-theory without separate 
valence quarks are trivially included in this. 

For the evaluation of the zero-mode integrals it is convenient to use the graded for- 
malism where partial quenching is achieved by adding additional bosonic and fermionic 
species. The zero-mode integral corresponding to the graded version of Eq. ( |3.2| ) with n 
bosons and m fermions is known in closed analytical form for an arbitrary mass matrix 

H, 

det[/i^ Jij+j-lilJ'i)]i,j=l,---n+m 



7^ 

n,m 



({/^i}) 



(A.l) 



where fXi = miTiV . Here ^'s are defined as J^u+j-iifJ-i) = (— for z = 1, • • • n 
and Ju+j-iifJ-i) = lu+j-iifJ-i) for i = n + 1, ■ ■ ■ n + m, where and Ki, are the modified 
Bessel functions. In this paper, we need the case with {n,m) = (l,iV/ -|- 1): 

^1,1+ Nf{^^b\^J'V, {fJ'sea}) 



i4) 



X det 



/ KM 

V 



luifJ-v) luifJ-l) IuifJ'2) ■■■\ 

fivlv+l{flv) fJ.2L+l{fJ-2) ■■■ 

fj-liu+2{fj"u) ^A^Iv+2{^J'l) A*2-^i^+2(/^2) ••• 



(A.2) 



Here ^f, = m;,Sy, /i^ = m^TiV and /ij = mi'EV, where nif,, and rrii denote the masses 
of the valence bosons, the valence quarks, and the physical sea quarks respectively. In the 
arguments, the set of sea flavors are denoted by {^isea} = {miSF, m2Sy, • • •}. 



When some of the quark masses are degenerate, one simply uses I'Hopital's rule. For 
example, for /xi = /X2 one has 

1 



^l_l+Ar^. (/^bl/^i;, {^sea})| 



Mi=M2 



/ Ki,{p.b) ivifj-v) Lip-i) jj-i^iu-iim) ■ ■ ■\ 

-fJ-bKiy + liflb) fJ-lJu + lifJ-v) fJ'lIu+lifJ'l) LifJ-l) 

f^bKi'+2ifJ'b) l^lli^+2il^v) fJ'lIu+2{fJ'l) fJ-lIu+lil^l) •• 



iNf 



X det 



V 



(A.3) 



7 



Partially quenched observables are computed by differentiating Eq. ( |A.2| ) with respect 
to suitable sources and subsequently taking the limit fib ^ f^v For example, the zero-mode 
integral of Eq. (|3.3|) is 



1 - d 

-{U^v + Ul^)u = ^l^ilJ'vAi^sea]) = " Hm Z^^^^^j^ ^{jJLb\lJiv, {l^sea}) ■ (A.4) 
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To obtain the result in a = (or even 6^0) QCD vacuum, one could numerically 
sum over topology with the weight given by the partition function above; 



or use the analytic expressions known for the Nj = 2, 3 cases |25]. 

B. Doing the zero-mode integrals first: a ^(1) toy model 
The perturbative scheme presented in this paper relies on the operator 

^[M\u-i)e+e{u^-i)M], 

being small, and more precisely of 0{p^) or smaller, for any value of the mass m as this mass 
is taken to zero at finite volume. The difficulty in giving a precise counting to this term lies 
of course in it being a combination of mass m, field .^(x) (both of which can be assigned 
clear countings) and the zero mode field U. It is therefore interesting that an alternative 
scheme exists, which gives identical results, but which assigns a definite magnitude to this 
term. This is achieved by first doing the zero mode integral exactly, and only subsequently 
performing the perturbative evaluations of the non-zero mode integrals^. In the general 
SU {Nf + A'^) case this becomes rather cumbersome in practice, but it has the advantage 
that all terms are explicitly ordered according to the expansion parameter 1/L. Here we 
illustrate it for the case of a simple U{1) toy model at the topological sector with = 0, 
where the zero-mode integrals are trivially performed. 
We thus consider the U{1) "chiral Lagrangian"^ 

£ = —d^U{x)d^U^x)--[m^e'^U{x)+h.c.] + ---, (B.l) 

where U{x) £ U{1). We separate into zero modes and non-zero modes, 

U{x) = [7e*v^«(")/^ = e^V2«o/^e^^«(")/^, (B.2) 

and do the analogue of fixing topology to = by integrating over 6. Expanding pertur- 
batively in the non-zero modes ^(x), we get 

£ = -^mcos{V2^o/F) (l - ^{xf/F^) + ^(a^e(^))' + • • • , (B.3) 



'^We are grateful to M. Liischer for suggesting this alternative method. 

*0f course, in the real Nf — 1 theory chiral symmetry is broken explicitly by the anomaly, and this 
Lagrangian is therefore not relevant for describing the the Nf = 1 theory. We use it only as a simple toy 
model to illustrate in a very transparent way the effect of first integrating over the zero mode. 
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We now perform the zero-mode integration over exactly with respect to the term shown 
to get the partition function 

Z = J mx)]Io (m^V {^-^j d^'xiixfy^ e-y^'-^^^^^' + ■■■, (B.4) 

up to overall irrelevant factors. 

Expanding the Bessel function using Iq{x) = Ii{x), and exponentiating the expanded 
term leads to the effective partition function 

Z = 7o(mSy) J [d^(x)]e- ^<^'^[i(^''«)'+^'^'«(^)'l + • • • , (B.5) 



where an effective volume-dependent pion mass is defined by 

7D — 

^ ~ F2 Io{mi:V)' 



(B.6) 



wc note that, as expected, this effective pion mass approaches the standard tree-level pion 
mass expression in the limit where mT,V oo: 

ml ImT.lF'^ = m\ as mST/ ^ oo, (B.7) 

while in the opposite limit we find 

277T,S 1 F'^V 

ml ^ -mT,V = ml as mEV 0, (B.8) 

Clearly, in the usual p-regime expansion this provides us with the standard massive 
propagator term, while in the chiral limit taken at finite volume V , there is no infrared 
problem due to the momentum sums being taken over non-zero modes only. When mY^V 
becomes of order unity we recover the usual e-regime expressions. 

But here we are interested in seeing how our expansion looks if we add and subtract 
the standard mass term, as is done in the main part of this paper. We therefore do a trivial 
rewriting, 

C = \{d,ax)f + lmU{xf 

= lid.ax)? + Iml^ixf + lsmU{x)\ (B.9) 



where 



6ml = ml -ml, (B.IO) 



Can we treat ^5ml^{xf' as a perturbation? Near the usual p-regime where mTiV —>■ oo 
we can use the asymptotic expansion of Bessel functions, 

to see that 

Sml = -^ + ---, (B.12) 
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This is of ©(p^) as expected in this U{1) theory. When instead mTjV is of order unity, we 
get 

6ml ml, (B.13) 

which is also of 0{p^). The point here is that we know the full analytical expression 

2mS 



Sml 



i?2 



/o(mSy) 



(B.14) 



for all values of m and V, and it is easy to check that this function is of 0{p^) everywhere. 
The term 5m^^(x)^ is thus explicitly found to be of NLO and we can treat it perturbatively. 

Of course, as already discussed in Section 5, away from the conventional p-regime 
expansion, separating out 6mlS,{x)'^ and treating it perturbatively amounts to a re-ordering 
and partial resummation of terms in this expansion. This is because the propagator is taken 
to be the conventional massive one even when ml is no longer of O^p"^), but smaller. The 
difference between a calculation based on this propagator and one where the propagator 
has been expanded in ml up to the needed order is always of yet higher order, and thus 
only illustrates the inherent uncertainty in any fixed-order perturbative expansion. 

C. Numerical evaluation of A(0,M2) and dM2A{0,M'^) 

The definition of 51 (M^) in Eq. (U) and 52(^2) = -9^/251 (M^), requires an infinite sum 
over the 4-vector = n^L^ where Li = L (i = 1, 2, 3) and L4 = T. In this appendix, we 
will suggest how to evaluate numerically gi(M^) and g2{M'^). 

For ML > 1, the expansion 



in terms of modified Bessel functions is useful, while for ML < 1, the polynomial expression 

167r"' ^ (n — 1)! 

n=l ^ ' 

using the shape coefficients /3nS Q is appropriate. Here 

o-f-^Yf 2 ^ a2-ln47r + 7-3/2 







yi/2j [v^H 



I * (w) « (w) - ' I > ' (C'^' 

S{x) = expi-Tik'^x), (C.5) 

fc=— 00 
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where 7 ~ 0.577215665 is Euler's constant and the summation in S{x) is typically well 
approximated by a truncation to |/c| < 20. 

Since the modified Bessel function Ku(x) is well-defined for the complex value of x, 
the both expressions above can be easily extended for the complex arguments using (m + 
\A^)^^^U.=u = {m^ + A2)i/4eiarctan(A/m)/2^ ^^le prescription given in Eq. gj). 

As shown in Fig. |8|, ignoring the contribution from |n^| > 20 in the summation in 
Eq. ( p.l| ), and n > 5 in the summation in Eq. ( |C.2| ), both formulas agree in a rather 
long interval around ML = 1. We also find an excellent agreement between these two 
expansions for both Kegi{M'^) and Re^2(-^^) at imaginary values of M^, using the same 
truncation. It is in any case trivial to increase the accuracy by including more terms. 

To summarize, we have used 



^ 5!(M2) ^ - M^ln(MVV2) _ Y^l^^ ^M2(«-i)y("-2)/2 (|M|L < 1) 



(C.6) 

and similarly for their derivatives 32- We plot Regi{MQ + iM^) and Reg2{MQ + iM^) for 
various choices of fixed Mq in Fig. |9|. 

Specifying the renormalization scale Hsub = 0.77 GeV of the low energy constants Lj's, 
and the above prescription for gi and 52 > 

A(0,m2) = ln^ + gi(M2), (C.7) 
9m2A(0,M2) = ln^ + 1 -92(M2), (c.8) 

V ^^ sub / 

can then be evaluated numerically. We note that in the small mass region the logarith- 
mic contribution InM^ is canceled by a similar one in gi and that both A(0,M2) and 
9^/2 A(0, M^) have no IR divergences in the limit — > 0. 

D. Useful properties of (^(a;, M^, M^) 

We also need the 1-loop correction from off-diagonals part of ^, G(0,M^,M^). It can, in 
principle, be expressed in terms of A(0,M^). In this appendix, we discuss the UV diver- 
gence of (5(0, M^, M^) and rewrite it using A(0, M^) and 9^/2 A(0, M^). Exphcit examples 
for the degenerate case and non-degenerate Nf = 2 + 1 fiavor theory will be given. 

First let us consider the UV-divergent part of A(0, M^), 
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Here the first term is quadratically divergent, the second term produces a logarithmic di- 
vergence, and the remaining 0{M*) terms are UV finite. As is well-known, the quadratic 
divergence is absent when we employ dimensional regularization. 



By expanding in the mass, the UV-divergent part of G(0,M^,M^) can be written 

1 



G{0,M\M^) = yYl 



1 ^ / 1 2M2 1 

Nf 

^ f f 



N 



-f \ 

^M^-^yMh\^ + ---, (D.2) 

where the logarithmic divergence of the last line is canceled by a renormalization of Lj's 
as seen in Section |3[ 

Although Eq. ( p^ ) shows that 0(0, M^,]^^) inherits the UV properties of A{0,M'^), 
the expansion is not useful when we want to see the finite part of (5(0, M^, M'^), since 
the omitted C(M'') terms are the at the same order as the former two terms. In order to 
obtain a convenient rewriting, let us define a function 



where k denotes the number of different quark masses and nj > 1 is the degeneracy of the 
i-th mass satisfying = Nf. Here we have ordered the masses Mf- < Mf_^^ -^-^ for any 

i. Noting f{t) is a monotonically decreasing function. 



and 

r2 j\/r2 



hm /(M2, + e) = oo, lim /(M2.i - e) = -oo, (D.5) 
fit) / 0, for t<Ml,Mlk<t, (D.6) 

one can show that an equation f{t) = has k — 1 different solutions (we denote them by 
t = Mfj), each of them satisfying 

Mf. <m2. <m2^i,+i, (l<i</c-l). (D.7) 
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We illustrate this in the plot Fig. 10 
Hence, 

-fi-p 



2, uH-Hp'+m 



and G(0, M^,M^) can thus alternatively be expressed as 



^f^t^o(p'+M')'ur\p'+m) 

fc-i 

J2 A(0, M^) + BA{0, M^) + C A(0, 



U'ip' + M; 



1 



where the coefficients Aj's, B and C are given by the residue of 

f2{t) 



(_t + M2)2nt-^(-t + M2)^ 



(Di 



(D.9) 



(D.IO) 



(or — (— t + M'^)f2{t) for C), at each pole. Note that C = when is equal to any of 
the physical masses. 



Noting that both A(0, M^) and A(0, M^) are infra-red finite in the limit — > 0, 



A(0,m2)|j;^2_o = 

9m2A(0,m2)|jv,2^o -- 

the chiral limit of G(0, M^,M^) is given by 

rfc-i 



Pi 



167r2 167r2 



(D.ll) 
(D.12) 



G(0,0,0) 



1 

w 



f 



E-4 



i|M2=0 



A(0,m2)_5| 



M2=0 



/3i 



(D.13) 



where UV divergence ci is absorbed into the renormalization of Lg- 

Here we give some useful examples. For the fully degenerate case, i.e., equal masses 
= M^g„ for all i, the above expression for G is greatly simplified. 



g(o,m2,m2) = ^ [A(0,m2) + (m2-MD5m2A(0,m2) 



(D.14) 



in agreement with the result presented in ref. |T^ ]. 

For an Nj = 2 + 1 flavor theory, the equation f{t) = is easily solved and one obtains 



G(0,M^m2) = - 



+ BA{0, M^) + C A(0, 



(D.15) 
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where M^^ = 2muE/F^ = 2mdS/F^, M^g = 2msS/F^ and the coefficients are given by 

_ - Mir ^^1^ 2{Mi-Mir 



^ ^ 3(M^-Mi)(M^-jV4) ^ 
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2 k 



partially quenched m--m^, \-0 
partially quenched m~m^, V-2 
unquenched m--m^^, V-0 
unquenched m— rh^, V-3 

'■■'»;!;— 



10 20 30 40 50 60 70 80 
m- ZV 

Figure 1: The function rriiTtV (^{Uu + Ujj}u/2 — ij is plotted in various cases in 2+1 flavor theory: 
A. partially quenched theory (m^ = my) with mu^d^y = 5,msT,V = 40 at — {rriu^d.s denotes 
the quark mass of u,d,s flavors respectively.), B. The same as A but at = 2, C. unquenched 
theory with m.^ — rriu.d at = where nisTiV — 40 is fixed. D. unquenched theory with = rus 
at v — 3 where rriu.d'^V = 5 is fixed. Every curve is kept at 0(1) at any value of m^, which confirms 
Ai{U — 1) always can be considered as NLO to the calculation of the chiral condensate. 



-27- 



Nj= 2, m„=10 MeV, L=2 fm 



A 

> 

1 > 
V 



• — , 
Q. 




new formula (v=0) 
p-expansion (0=0) 
£-expansion (v=0) 



Nf = 2, L=2 fm, m^=10 MeV 




Figure 2: The chiral condensate (top) and Dirac spectral density (bottom) for the case with 
m„ — rud — 10 MeV, L — T 12 = 2 fm. The result of this paper is given by the thick line which 
smoothly connects the one in the e-regime (solid) with the on of p-regime (dotted). Parameters 
have been chosen as E^/^ ^ 250 MeV, F = 90 MeV, Lg(/i^„b = 0.77GcV) = 0.05 x lO^^, Lg(^,„b = 
0.77GeV) = 0.5 x lO^^ and HUpLsub = 0.77GeV) = 0.1 x lO-^. 
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Figure 3: The same plots as in Fig. the condensate (left) and the spectral density (right) 
but with the larger volume size L — T /2 — i fm. The convergence of three expansions around 
MvvL ~ 2-3 or XT,V ~ 10 becomes better. 




Figure 4: The condensate (left) and the spectral density (right) at fixed topology v — for 
different sea quark masses at 2 MeV (solid), 10 MeV (dotted), 30 MeV (small dotted) are given. 
We again set L ^ T/2 = 2fm. The black filled circle in the left panel shows the physical points 
where my = niu — riid- 
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Figure 5: The topology v dependence of the condensate (left) and the spectral density (right) at 
a fixed value of to„ = 10 MeV and L = T/2= 2fni and H^{Q.nGe\f) = 0.1 x 10"^^ The topology 
dependence becomes negligible around Aly^L = 3-4. 




Figure 6: The (full theory) condensate (left) where — to„ = irid and the spectral density (right) 
at fixed m„ = irid = 10 MeV, are given at 6* = 0. Due to the ambiguity of H2, the condensate is not 
a well-defined quantity when iW^ or Myy ^ 0, while the spectral density is free from this ambiguity. 



-30- 



m, =10 MeV, L=2 fm 





Figure 7: The condensate (top) and the Dnac spectral density (bottom) with Nf = 2 + 1 theory 
are given (soUd curve). For the comparison, Nf — 2 plot is also presented (dotted). As the valence 
mass or the eigenvalue increases, both of them show larger deviation from 1 than the case with 
Nf = 2. 
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Figure 8: The plot of gi{M'^)V'^/'^ (left part) and g2{M^) (right part) for T = 2L. The truncated 
sums in Eq. ( |C.6| ) for 2i^'^) (solid line) and §{ 2(Af^) (dotted) agree well around ML ^ 1. Note 
that 5i(M2)1/1/2 converges to -/9i = -0.08360 in the chiral hmit. Note also that |.gi(M^)| attains 
its maximum at M = 0. 




Figure 9: Plot of Regi{M§ + iM^)V^/^ (left part) and Reg2{M§ + iJVP) (right part) for T = 2L. 
for three fixed values of MqL = (solid), 0.05 (dotted) and 0.1 (small dotted). 



-32- 



-0.5 0.5 1 1.5 

t 

Figure 10: Example of f(t) in Eq. ( |D.3| ), where we have N f = (!>, k — b, rii = l(z ^ 2), n2 = 2, 
and {Mf^\ = {0.1,0.3,0.5,0.75,0.9}. One can see fit) = has fc - 1 = 4 solutions between the 
poles (plotted by the crosses). 
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